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I present a method for evaluating the hadronic vacuum polarization contribution below 1 GeV to a^j which relies 
on analyticity, unitarity and chiral symmetry, as well as on data. The main advantage is that in the region just 
above threshold, where data are either scarce or have large errors, these theoretical constraints are particularly 
strong, and therefore allow us to reduce the uncertainties with respect to a purely data-based evaluation. Some 
preliminary numerical results are presented as illustration of the method. 



1. Introduction 

The anomalous magnetic moment of the muon 
is now known experimentally to an extremely 
high precision , and can therefore be used as a 
thorough test of the standard model. The most 
uncertain and debated part of the standard model 
calculation of is the hadronic vacuum polariza- 
tion contribution aJl'^P . In order to make this test 
of the standard model a significant one, we need 
to make substantial progress in the evaluation of 
the hadronic contributions, and in particular of 
the leading one, - which was precisely the 
scope of this Workshop. This contribution can 
be calculated in terms of another experimentally 
measured quantity, the cross section e~^e~ — s- 
hadrons, and any improvement in the measure- 
ment of this cross section is immediately reflected 
in aJj^P . Indeed many discussions at the Work- 
shop concerned different possible ways to measure 
the e~^e~ hadronic cross section, and have given 
us an overview of an impressive amount of ex- 
perimental and theoretical work devoted to this 
problem. Given a set of data points at different 
center of mass energies for a{e^e'' — > hadrons), 
the evaluation of ai^^P amounts to the calculation 
of an integral of a function once one knows it at a 
discrete set of points. This problem is standard, 
and can be solved in different ways - the simplest 
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of them being the use of the trapezoidal rule, as 
done, e.g. in Refs. i2|3,. 

The use of the trapezoidal rule, or variants 
thereof, reduces the role of theory to an abso- 
lute minimum ~ the evaluation of a'^^P is done 
almost only with data. The role of theory can- 
not be reduced to zero because it is needed in the 
two extreme regions s ~ 4M^ and s — s- cx). In the 
latter region one can make use of perturbative 
QCD, and in the former one of chiral perturba- 
tion theory (CHPT). Actually, CHPT cannot give 
a sharp numerical prediction for the pion vector 
form factor Fv{s) (the dominating contribution 
below 1 GeV) close to threshold: it only predicts 
that the behaviour is very smooth, well approx- 
imated by a polynomial, whose coefficients are 
related to a few of the low energy constants of 
the chiral Lagrangian. The use of CHPT in this 
region amounts to an extrapolation of the data at 
somewhat higher energies down to threshold with 
a polynomial of low degree 2 3 . 

My aim here is to show that not only close to 
threshold, but even up to 1 GeV the use of some 
theory is actually quite useful, especially because 
it allows one to reduce the uncertainty in - 
and this region contributes the largest fraction of 
the total error. Theory in this case means some 
very general properties which we know the vector 
form factor -FV(s) must satisfy: analyticity and 
unitarity. Gombining these properties with chiral 
symmetry, which is relevant in the very low en- 
ergy region, one can construct a representation 
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which is very constraining for the form factor: 
the remaining httle freedom can be fixed with the 
help of data as I illustrate in what follows. 

2. Definition of a^^P 

The leading hadronic contribution to {g — 2)^^ 
is due to the hadronic vacuum polarization cor- 
rection shown in Fig. ^ 




Figure 1. Hadronic vacuum polarization contri- 
bution to (g — 2)^. 



This contribution to {g — 2)^ is of order and 
can be expressed in terms of the cross section 
crfe+e" 



a H: 



hadrons) evaluated to leading order 



V 37r y 



as 



/4Af2 S 

where K{s) is a known kernel, see e.g. and 
r(o)(g+g- ^ hadrons) 



(1) 



47ra2/3s ' 

where the superscript indicates that the cross sec- 
tion has to be taken to leading order in a^. At 
low energy ^/s < 1 GeV, the contribution of the 
two-pion state dominates the cross section 
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R('\s) ^ R-^s) ^ i (1 - 4M^/sf' \Fy{s)\' .(3) 

^We do not discuss here the problem of extracting the 
leading order term from data. 



The latter, as shown, is given by the vector form 
factor of the pion (again to leading order in a), 
which is the quantity we will now discuss. 

3. The pion vector form factor 

The pion vector form factor ^V('S) is an ana- 
lytic function of s in the whole complex plane, 
with the exception of a cut on the real axis for 
s > 4M^: approaching the real axis from above 
the form factor stays complex and can be de- 
scribed in terms of two real functions, its modulus 
and phase: 



(4) 



Omnes ^ has shown that analyticity relates the 
modulus and the phase, such that the whole func- 
tion can be given in closed form in terms of its 
phase: 



Fvis) = P(s)exp 



dx ■ 



Six) 
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(5) 



where P(s) is a polynomial which determines the 
behaviour of the function at infinity, or number 
and position of its zeros. In order to respect 
charge conservation -FV(O) = 1, we set P(0) = 1. 
The representation ^ makes it apparent that if 
one knows the phase on the cut and the zeros of 
the form factor one can calculate the form factor 
everywhere in the complex plane. 

In the elastic region Watson's theorem relates 
the phase of the vector form factor to the phase of 
the TTTT scattering amplitude with the same quan- 
tum numbers, I = i = 1: 



5{s) — Sl{s) for s < Si- 



16M' 



(6) 



While the inelastic threshold y^Sm — .56 GeV ap- 
pears to be rather low, inelastic contributions in 
this channel are known to become relevant only 
above the KK threshold. To an excellent ap- 
proximation the phase of the vector form factor 
coincides with the tttt phase shift Si up to 1 GeV. 
The latter has been recently studied in the frame- 
work of Roy equations, first with no extra input 
|7], and also in combination with chiral symme- 
try jH). The upshot of these analyses is that Si 
is constrained to a remarkable degree of accuracy 
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up to about 0.8 GeV. In our approach we make 
explicit use of these general properties of the vec- 
tor form factor, and of our knowledge of the tttt 
phase shifts. Our strategy can be summarized in 
the following points: 

1. we construct a convenient representation 
which automatically respects the properties 
of analyticity, unitarity and chiral symme- 
try; 

2. we fix the free parameters which appear in 
this representation by fitting data; 

3. we evaluate the integral CQ) up to 1 GeV us- 
ing our analytic representation of the form 
factor. 

A similar strategy was also adopted in Ref. [5], 
but with a limited use of chiral symmetry. 

4. A convenient representation of the vec- 
tor form factor 

As discussed above, the vector form factor in 
the low energy region is to a large extent dom- 
inated by the tttt phase shift 61- Inelastic ef- 
fects, which are small, but nonnegligible at the 
needed level of accuracy, will also be taken into 
account and parametrized in terms of a smooth 
function which has the correct analytic proper- 
ties. In order to evaluate a^^^ we need the vector 
form factor to leading order in a, i.e. with elec- 
tromagnetic interactions switched off - our form 
factor Fy does not include vacuum polarization 
corrections. To make the connection to the tttt 
phase determined in [Hj it is also convenient to 
work in the isospin limit of strong interactions^ 
rriu — rud- In the tttt scattering amplitude these 
isospin-violating effects only show up at order 
(m„ — m^)^ and are negligible. In the form factor, 
however, they are linear in the quark mass differ- 
ence, and moreover they are enhanced, in a cer- 
tain energy region, by the small mass difference 
between the p and w mesons, which appears in the 
denominator. This enhanced isospin-violating ef- 
fect cannot be neglected at the level of accuracy 

^The quark mass should be chosen such that the common 
pion mass is equal to the physical charged pion mass 



which we are working at. We are therefore going 
to represent the form factor as a product of three 
functions that account for the prominent singu- 
larities in the low energy region (10): 

F{s) ^ Gi{s) ■ G2{s) ■ GUs) ■ (7) 

The first term is the Omnes factor that de- 
scribes the cut due to 2tt intermediate states: 

f s dx5\{x) \ 

Gi{s)=expi- — -\ . (8) 

[ 71" Jam^ x[x- s) J 

The phase Si which enters Gi (s) is obtained with 
an updated version of the analysis in |Hj : in solv- 
ing Roy equations we need an input for the imag- 
inary part of the various partial waves at and 
above 0.8 GeV. In [S] the input I = i = 1 partial 
wave had been fixed with the CLEO data JI] on 
the vector form factor: the most important input 
parameter, the phase at ^/sq = 0.8 GeV had been 
taken equal to 

Ji(so) = (108.9±2)°. (9) 

Above that energy we had used the parametriza- 
tion of Hyams et al. ^21- In the present work 
we leave (^i(*o) and Sl{si), with y/si = 1.15 GeV 
(the upper limit of validity of the Roy equations), 
as free parameters and determine them by fitting 
data on the form factor. Once these input param- 
eters are given, Roy equations and chiral sym- 
metry fix uniquely the phase between these two 
points and all the way down to threshold, as il- 
lustrated in Fig. [3 The phase shown in the plot 
corresponds to 

,51 (so) = 109.5° Slisi) = 165.9° (10) 

which are typical values we obtain in our fits, and 
in perfect agreement with the input phase used in 
Uj. The error which we get on (5x(so) is however 
about a factor two smaller than that in 

The function Gi^(s) contains the pole generated 
by Lo exchange, 

G^(s) = 1 + e — ^ + --- 

= (M^-iir„)2. (11) 

The pole term cannot stand by itself because it 
fails to be real in the spacelike region. We re- 
place it by a dispersion integral with the proper 
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Figure 2. The tttt phase shift 6l: the value of 
the phase at the two (red and blue) dots are free 
and the Roy equations and chiral symmetry com- 
pletely fix the phase everywhere else. 



behaviour at threshold, but this is inessential: the 
representation for G^{s) that we are using is (a) 
fully determined by the values of e, M^^ and 
and (b) in the experimental range, |Gr^(5)| is nu- 
merically very close to the magnitude of the pole 
approximation. 

The function G2(s) represents the smooth 
background that contains the curvature gener- 
ated by the singularities not acconted for by Gi 
and Goj. We analyze this term by means of a 
conformal mapping. The Air channel opens at 
Sin = 16 M^, but phase space strongly suppresses 
the strength of the corresponding branch point 
singularity, which is of the form (1 — Si„/s)^/^. 
The transformation 



y Sin — Si + \/ Sin ~ S 

maps the plane cut along s > Sin onto the unit 
disk in the ^-plane. It contains a free parameter 
si - the value of s that gets mapped into the ori- 
gin. We find that if si is taken negative and suffi- 
ciently far from the origin, the fit becomes insen- 
sitive to its specific value. We set si = — 1 GeV^. 
We approximate G2(s) by a polynomial in the 
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Figure 3. Typical shapes of the factors G2 and 
Guj that enter the vector form factor. 



variable z, 

p 

G2{s)^l+J2c^{z'-4), (13) 

4=1 

where Zq is the image of the point s — 0. The 
terms involving zq ensure that the background 
does not modify the charge, 62(0) = 1. The con- 
dition that the branch point singularity has the 
form (1 — Si„/s)^/^ implies four constraints on the 
coefficients: if we want a nontrivial contribution 
from G2(s) we need at least a fifth-order polyno- 
mial. In the following we will vary P, the order 
of the polynomial, between and 8. Both fac- 
tors G2 and Guj are very small in the region up 
to 1 GeV, as illustrated in Fig. 0| only G^j gener- 
ates a 10% effect but only in a very narrow region 
around M^^, otherwise both give an effect at the 
few percent level. Notice that in principle G2 may 
have zeros somewhere in the complex s plane, as 
allowed by the general Omnes representation 
- we make no assumptions on their position, nor 
on their number, if we vary P, and let the data 
choose where they should lie. 

5. Numerical results 

Our representation of the vector form factor in- 
volves three functions, each of which has free pa- 
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xVd.o.f. 


XcMD2 


XnA7 


IQiOflp 









84.9/83 


43.6 


43.7 


420.1 ±2.1 


489.5 ±2.2 


0.4254 ±0.0020 


5 


78.4/82 


35.9 


42.6 


423.8 ±2.6 


494.1 ±2.7 


0.4300 ±0.0024 


6 


78.1/81 


36.0 


42.2 


424.4 ±2.8 


494.7 ±2.9 


0.4339 ±0.0051 


7 


73.5/80 


31.7 


42.2 


423.4 ±2.9 


493.2 ±3.0 


0.4350 ±0.0051 


8 


73.5/79 


31.6 


42.2 


423.5 ±5.7 


493.4 ±7.4 


0.4347 ±0.0052 



Table 1 

Numerical results for fits to CMD-2 [TSI and (spacelike) NA7 data [HI- The CMD-2 data used here do 
not contain vacuum polarization effects nor final state radiation. The errors given are purely statistical. 



rameters which we pin down by fitting data. The 
free parameters at our disposal are the following: 

• in Gi we have two free parameters, the 
value of the phase 5\ at 0.8 and at 1.15 GeV; 

• in G2 we have P — A free parameters, where 
P is the degree of the polynomial in the 
conformal variable z. We will vary P within 
a reasonable range of values, and see how 
the results depend on it; 

• in Guj we have in principle three free pa- 
rameters, Ei^, and ~ however, since 
the mass and width of the to are rather well 
known from other experiments we fix them 
at the PDG values. We allow, however, the 
energy calibration (to which G^j is very sen- 
sitive) to shift within the estimated experi- 
mental systematic uncertainty. 

All in all we have 4 + (P — 4) free parameters, 
depending on the degree of the polynomial used 
to describe inelastic effects. 

An example of our numerical results obtained 
by fitting the most recent CMD-2 data and the 
spacelike NA7 data ^31 for values of P between 
and 8 is given in Tabled The quantities de- 
noted by Op and a2A/jf correspond to the integral 
(0 cut off at 0.81 GeV and 2Mk respectively, 
whereas (r^) is the square of the pion charge ra- 
dius. It is worth stressing that the good fit ob- 
tained with P = shows how well the two-pion 
intermediate state alone describes the behaviour 
of the form factor both in the timelike and the 
spacelike region. The numerical results obtained 
in the first line can be considered as a theoretical 



prediction of Up based on the calculation of the 
TTTT phase shift done in jS] (of course the results 
in that paper is also based on some experimen- 
tal input). It is clear, however, that at the level 
of precision needed for this prediction is not 
sufficiently precise and we must explicitly account 
for the inelastic effects encoded in G2. 

If we switch on the function G2 and allow for 
one free parameter in it (P = 5) we observe a size- 
able improvement in the x^, mainly in the part 
which comes from CMD-2 data - the addition of 
a further parameter (P = 6) slightly improves 
the fit to the NA7 data. With P = 7 we have 
again a sizeable decrease in Xcmd2j whereas go- 
ing to P = 8 brings no improvement anymore 
but leads to a substantial increase in the errors 
on all calculated quantities. Going to even higher 
values of P does not make sense anymore, and 
it is reasonable to choose as final result the one 
with P — 7, while the variation of the result with 
P will be included in the final uncertainty, e.g. 
by taking the difference among the P = 6 and 
P — 7 results as theoretical uncertainty. These 
numbers are preliminary and given only for illus- 
tration purposes: the most important point to 
stress concerns the uncertainties, rather than the 
central values. A comparison to results obtained 
with the trapezoidal rule, like e.g. the recent up- 
date of Jegerlehner'* 5 

Up = 429.02 ± 4.95 (stat.) (14) 

shows that the reduction in the purely statistical 
error is substantial. Final results will be given in 
a forthcoming publication 15. We plan to extend 

*Tho difTeronce in the central value is mainly due to final 
state radiation, included in 1141 but not in Tab. [Tl 
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Eniax/KnaxiGeV) 


XcMD2 


XoLYA 


Xaleph 


XcLEO 


421.5 




.81/- 


20.3(27) 


14.1(26) 






431.0 




-/.81 






8.1(10) 


16.3(21) 


427.2 




.81/.81 


21.1(27) 


22.4(25) 


11.0(10) 


19.7(21) 


427.4 


497.7 


.97/.81 


38.2(43) 


27.3(42) 


12.0(10) 


20.1(21) 


427.7 


501.0 


.81/.97 


22.0(27) 


23.9(25) 


12.3(16) 


31.3(28) 



Table 2 



Results of simultaneous fit to e~^e~ and r data in different energy regions. The third column gives the 
maximal energy up to which the e"'"e~ or r data are fitted. The first two rows concern fits done to cither 
e'^e~ or r data. The number in brackets next to the value gives the number of data points fitted. 



the analysis to other sets of data, like older e+e~ 
data (e.g. jJH), but also data on the weak vector 
form factor from r decays |17lll| . 

As is well known, the latter sets of data show a 
systematic deviation from the e'^e~ ones, even af- 
ter correction for known isospin violating effects 
|18| . One should stress that the calculation of 
these effects, although done on a sound theoret- 
ical basis, cannot be improved systematically: it 
is difhcult to estimate the final theoretical uncer- 
tainty in the calculation of these corrections, and 
it is therefore not possible to exclude that the re- 
maining discrepancy between e+e~ and t data 
is actually due to an unaccounted isospin violat- 
ing effect. Indeed this possibihty (in particular 
a difference in mass and width between charged 
and neutral p mesons) has been discussed at this 
Workshop and in a recent publication [501 . 

One of the striking features of this discrep- 
ancy is that it has a peculiar energy dependence: 
as observed in [2], the two sets of data are in 
good agreement below about 0.8 GeV, but show 
a marked difference from that energy on. It is 
in fact possible to make a good fit to both sets 
of data if one limits oneself to the energy region 
below 0.8 GeV, as illustrated in Tab.El The first 
two rows give the results of the fits to either the 
e+e~ or the r data: the outcome for Op shows a 
discrepancy of about 10 units, which is certainly 
larger than the error one would like to achieve in 
this determination. In the third row, however, it 
is shown that if one fits simultaneously all data 
sets one can still obtain a good for all data 
sets, and a value for Op which sits between the 
e'^e~ and r value, somewhat closer to the latter. 



The last two rows show that if one extends the 
fit region higher up only for one of the two sets 
of data, the result for Op remains stable, whereas 
the discrepancy in the value of a2M^ is reduced 
to about three units. 



6. Conclusions 

I have discussed a method of calculation of the 
low energy contribution to which relies as 
much as possible on theory - in the form of an- 
alyticity, unitarity and chiral symmetry. These 
properties do constrain the vector form factor of 
the pion below 1 GeV quite strongly and can be 
used to safely interpolate between data. In par- 
ticular, it is sufficient to have very precise data 
in the p region to pin down the free parameters 
in the form factor and make a controlled extrap- 
olation down to threshold. In this manner one 
can sizeably reduce the error in aJ^/P coming from 
this energy region, where data are scarce and with 
larger errors. I have illustrated this with a few nu- 
merical examples and stressed, in particular, that 
because of the agreement between e^e~ and t 
data below 0.8 GeV, one can give a stable predic- 
tion for the region below this energy. The analysis 
is still in progress, and final results will be given 
in a forthcoming publication |15j . 

In this Workshop we have heard of plans or 
ongoing efforts for measuring the pion vector form 
factor in several different ways, which means that 
in a few years time there will be several new sets 
of data to which this method can be applied. It 
will be extremely interesting to look at the picture 
that will emerge from these. 
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